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ABSTRACT 

In this paper we present a new method for dating the surface of the Moon, obtained by 
modehng the incoming flux of impactors and converting it into a size distribution of resulting 
craters. We compare the results from this model with the standard chronology for the Moon 
showing their similarities and discrepancies. In particular, wc find indications of a non-constant 
impactor flux in the last 500 Myr and also discuss the implications of our findings for the Late 
Heavy Bombardment hypothesis. We also show the potential of our model for accurate dating of 
other inner Solar System bodies, by applying it to Mercury. 

Subject headings: solar system: general — planets and satellites: Earth, Mercury, Moon 



1. Introduction 

Craters are among the most spectacular surface 
features of the solid bodies of the Solar System. 
Cratering studies provide a fundamental tool for 
the age determination of planetary and asteroidal 
terrains. Since the beginning of the lunar explo- 
ration, age estimates for the lunar terrains were 
derived, followed by chronology models for the 
other terrestrial planets. The development of the 
lunar chronology greatly helped in interpreting 



the evolution of the Solar System and in particu- 
lar of our own planet, the Earth. Recently, thanks 
to a fleet of new space missions (Mars Express 
to Mars, MESSENGER to Mercury, and Kaguya 
to the Moon, only to name a few), this field of 
research entered a new exciting phase, where ac- 
curate age estimates provide means for detailed 
small-scales geological studies. 
The method currently used for dating purposes 
defines a chronology (crater surface density vs. 
absolute age) for a reference body for which ra- 
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diometric ages are available for different terrains. 
Then, on the basis of models predicting the im- 
pactor flux ratio between the reference body and 
another generic body, it is possible to estimate 
the age of the latter. This method can be defined 
experimental, since it develops a chronology of a 
reference body for which two measurable quanti- 
ties are available: absolute ages and number of 
craters for selected areas. The Moon is the only 
body of the Solar System that could fulfill both 
requirements so far. This method has been de- 
veloped and refined in the last 30 years by many 
researchers and r epresents the reference for dat- 
ing p urposes (e.g. Hartmann et al. 19811 iNeukum 
19831 ). In this approach, the development of the 
chronology for the reference body does not depend 
explicitly on the physics of the cratering process. 
However, such information becomes necessary -via 
a cratering scaling law- in order to apply the Moon 
chronology to other bodies or to infer the flux of 
impactors from the observed cratering record. 
A possible alternative -purely theoretical- ap- 
proach could be based on the accurate estimation 
of the time-dependent impactor flux for each tar- 
get body, through the history of the Solar System. 
With this method, a direct comparison between 
the observed distribution of craters on a given ter- 
rain and the isochrones produced by the model 
would give its age. Although the knowledge of 
the formation and evolution of the Solar System 
greatly improved in recent years, such ambitious 
goal is presently far beyond the capabilities of the 
available models, due to the large uncertainties in 
the early stages of its formation. 
In this paper we introduce a third approach for 
determining the chronology of objects in the inner 
Solar System. This method, which can be con- 
sidered as a hybrid of the first tw o approaches 



is bas e d on t he dyn amical model bv lBottke et al 



(|200d . I2OO2I l2005al |bh which describes the forma- 
tion and evolution of asteroids in the inner Solar 
System. In this framework, first the flux of im- 
pactors on the Moon is derived, then the impactor 
size distribution is converted into a cumulative 
crater distribution via an appropriate scaling law, 
and finally the radiometric ages of different regions 
-the Apollo and Luna landing sites- are used for 
the calibration of the lunar chronology. The main 
advantage of our approach is that the adopted 
model naturally describes the distributions of the 



impactor sizes and velocities for any body in the 
inner Solar System. Thus, the lunar calibration 
can be exported with precision to any other body. 
A drawback of the method, however, is that the 
chronology of the reference body depends -through 
the scaling law- on the physics of the cratering 
process, which is as yet not fully understood. Fur- 
thermore, the precision of the age estimate relies 
on the accuracy of the dynamical model. Encour- 
agingly, the adopted dynamical model is capable 
of providing a good representation of the present 
asteroid population and size distribution both in 
the near-Earth space and in the main belt, with 
a maximur n deviation a.mounti ng to less than a 



factor of 2 (jBottke et al.ll2005bl ) 



2. Modeling the cratering on the Earth- 
Moon system 

The present inner Solar System is continu- 
ously reached by a flux of small bodies having 
sizes smaller than few tens of km. A fraction 
of these bodies (namely those having perihelion 
< 1.3 AU) are called near-Earth objects (NEOs). 
Contributors to this flux are represented by main 
belt asteroids (MBAs) and Jupiter family comets 
(JFCs). Both contributions have been modeled 
by a number of authors and it is generally ac- 
cepted that the MBAs constitute the main source 
for the flux presently observed in t he inner Solar 
System (e.g. iMorbidelli et al.l 120021 ) . therefore we 
neglect the cometary contribution. This flux is 
sustained by a few fast escape tracks, which are 
continuously replenished with new bodies as a re- 
sult of coUisional processes and semimajor axis 
mobility. The main gateways to the inner Solar 
System are the vq secular resonance with Saturn 
and the 3:1 mean motio n reson a nce with Jupiter 



Morbidelli et al 



(IMorbidelli k GladmanI Il998 
2002) • One of the most recent and accurate models 



co ncerning NEOs formation has be en developed 
bvlBottke et all (|2000l . l2002h . while iBottke et af 



(|2005aM bl) modeled the main belt dynamical evo- 
lution. These models have been adopted in order 
to estimate the properties of the impactor flux at 
the Moon. 
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2.1. Modeling the impactor flux: the sta- 
tionary case 

The flux of impactors may be described as a 
differential distribution, (f)(d,v), which represents 
the number of incoming bodies per unit of im- 
pactor size (d) and impact velocity (v). An im- 
portant quantity -since it can be constrained from 
observations- is the size distribution of the incom- 
ing bodies, namely the number of impactors per 
unit impactor size. Let h{d) be such distribution, 
we then have: h{d) = J <j){d, v) dv. Therefore, 
without loss of generality, we can write: 



cl){d,v)^h{d)f{d,v) 



(1) 



where f{d,v) is the distribution of the impact 
velocities (i.e. the impact probability per unit im- 
pact velocity) for any dimension d, and which ful- 
fills the normalization constraint / f{d,v) dv = 1. 
We used the NEO orbi t al distributi on model com- 
puted bv lBottke et~an (|2000l [20o3 ). generated by 
integrating two sets of test particles placed into the 
Vq and 3:1 resonances (with 2 x 10^ and 7x 10^ par- 
ticles, respectively). For each particle they com- 
puted the impact probability and the impact ve- 
locity r elative to the targe t body. Following the 
work o f lMarchi et £ll(|20ok ). we computed the im- 
pact velocity distribution (corrected for the gravi- 
tational cross-section of the targets) . We assumed 
that the Moon revolves around the Sun with an 
orbit identical to that of the Earth. Moreover, 
we corrected the impact probability provided by 
the Bottke et al model in order to account for 
the different collisional life-time of the particles. 
This is necessary because the orbital evolution of 
the particles is coupled wi th their size through 



the relationship r oc Vd ( Farinella et al. 19981 : 



Bottke et alll2005bl ). with T being the particles 



collisional life-time. In figure [T] we show the re- 
sulting impactor velocity distribution f{d,v) for 
the two extreme size bins used in our simulations 
and for both the Moon and the Eartt0. 



^The overall shape of the distributions vaguely resembles a 
Maxwellian distribution. It is interesting that the distribu- 
tions have a large number of high velocity (i.e. > 20 km/s) 
impactors. This distribution may help to explain some 
anomalies found for some terrestrial impact craters. In 
particular, only 3 out of about 60 known craters with ages 
< 60 Myr have been at present associated to tektite strewn 
fields l|French I Il998l ) . Since tektites are highly shocked 
super-heated melts ejected during hypervclocity impacts 



The impactor size distribution for the Earth (/ig) 
can be written as follows: 



he{d) = Pehnid) 



(2) 



where /i„ is the NEO differential size distri- 
bution and Pp is the intrin s ic colli sion probabil- 
ity Following iBottke et al.l (|2005bh . we adopted 
Pe = 2.8 • 10~^ yr~^. The distribution /i„ has been 



derived from iBottke et al.l (|2002i ). and it can be 
scaled to the Moon and other bodies. In the case 
of the Moon, this is done considering the Earth- 
Moon ratio of th e impact probability as a function 



Moon ratio or tn e impact proDaDuity as a tunction 
of impactor size (iMarchi et al.ll2005h . Let S(d) be 



this ratio, then the Moon impactor size distribu- 
tion becomes: 



hmid) = PeK{d)Y.{d) 



(3) 



In figure [2] the cumulative impactor size dis- 
tribution for both the Moon and the Earth are 
shown. A similar approach can also be used for 
evaluating the fiux on other targets, like Mercury 
(see section [6]). In equation [1] the dependence of h 
and / over time is neglected since we are consider- 
ing a stationary problem, while eq. [2]and[3]strictly 
describe the present flux. Arguments about possi- 
ble variations in the flux intensity and shape have 
been discussed since the late 70s, and this still re- 
mains one of the most debated topics in the field 
of lunar crater ing. In t heir work , Neukum and 
coworkers (e.g. lNeukum| [l983: Neukum fc Ivanov 
19941 ) generally assume a constant shape of the 
impactor size distribution over time, while others 
claim a sudden change in the shape after the Late 
Hea vy Bombardment ( LHB), about 3.8 Gyr ago 
(e.g. Strom et al. 20051 ). One of the potentials of 
our model is that we can implement and evaluate 
the effects of a time-dependent size-distribution 
on the crater cumulative distribution, as shown in 
section [5] 

2.2. The scaling law 

In order to convert the flux of impactors into 
a resulting crater distribution, wc need to apply 



(e.g. lGlasslll990l : iKoeberl et alJll996^ . from simple statis- 
tical arguments we may estimate that a threshold perpen- 
dicular impact velocity of at least 24 km/s is required for 
their formation (more than 33 km/s modulus for a 7r/4 im- 
pacting angle). 
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a so-called crater-scaling law. Such scaling law 
attempts to describe the outcome of a cratering 
event based on the impact energy and on the phys- 
ical properties of both the target and the projec- 
tile. Once an appropriate scaling law is chosen, we 
may write the final crater diameter D as a function 
of the impactor diameter, impact velocity, and a 
number of parameters (density, strength... indi- 
cated by p), namely D — S{d,v,p). Despite the 
great effort (both computational and experimen- 
tal) devoted to this task, the physics of crater for- 
mation is far from being completely understood. 
For this reason, we decided to explore how the pre- 
dicted crater distribution depends on the choice of 
the scaling law. For this purpose we considered the 
three sc aling laws most used i n this fi eld. The first 
one by iHolsapple fc Housen ^M) (H&H here- 
inafter), reads: 



D^kd 



2vl 



Y 



(4) 



where g is the target gravitational acceleration, 
is the perpendicular component of the impactor 
velocity, 5 is the projectile density, p and Y are 
the density and tensile strength of the target. The 
quantities k and /i depend on the target material 
and V on its porosity. H&H estimated these latter 
parameters by best-fitting over a range of exper- 
iments done with different materials. From this 
study we adopted the values k = 1.03 and /i — 0.41 
for cohesive soils, and k = 0.93, p — 0.55 for rocks. 
We used v = 0.4 in all cases (H&H). Equation H 
accounts for the transition between strength and 
gravity regime, allowing a smooth transition be- 
tween those extreme conditions. 
The second scaling law considered is re ported by 
Ivanov et al.l ( 200lh . (adapted from .Schmidt fc Housen 



19871) (I&S hereinafter) 



D 



1.28 



d{5/p) 



0.26, ,0.55 



[{Ds, + D)gf 



28 



(5) 



where D^g represents the strength-gravity tran- 
sition crater. Dgg has been set equal to 120 m and 
30 m for the Moon and the Earth, respectively 
1996h . 

have co nsidered the scaling 
2004i) 



(jAsphaug et al 



Finally, 
used _b^ 



.Stuart & Binz el , - - ^ 

Shoemaker et al.lll990f ) (S&S hereinafter) 



law 

(adapted from 



D = 0.0224 W 



0.294 
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1/6 



(sina)2/3 (6) 



where, W is the impactor kinetic energy, is 
the Earth gravity and a is the impactor angle with 
respect to the surface (vertical impacts correspond 
to a = 7r/2). The numerical multiplicative factor 
takes into account the c orrection from tr ansient 
to final crater dimension ( Stuart fc Binzell 2004) . 
The crater size D reported in these formulas 
should be regarded as th e size of the final crater. 
However, following |Pi ki (Il980l) . a correction for 
the transition between simple and complex craters 
has also been applied to all the scaling laws in the 
form of: 



D = 
D = 



D 



if D<D^ (7) 
if D>D, (8) 



where Z)* is the diameter of the transition from 
simple to complex crateJl. All impacts are as- 
sumed to occur at the most probable impact an- 
gle, namely 7r/4 with respect to the normal to the 
surface. 

All the scaling laws assume that the density and 
the strength of the target are constant through- 
out the body. However, most target bodies are 
characterized by a layered structure, with a den- 
sity increasing with depth. Therefore, in order to 
produce a better representation of the target den- 
sity, we have computed the average density over a 
depth of about 10 times the radius of the projec- 
tile. This is because the size of the crater is -on 
average over different impact conditions- about a 
factor ten larger than the size of the impactor. 
In particular, in the case of the Moon we have as- 
sumed a 10- km layer of fractured silicates (megare- 
golith and heavily fractured anorthosites), on top 
of a bulk anorthositic crust in turn laying above a 
peridotitic mantle (see fig. [3]). For the Earth (see 
fig- [3|) w e considered the following lithospheric 
structure ( Andersonl l2007t ) : a 2 km thick layer 
of sedimentary rocks, a mainly granitoid upper- 
crust down to a depth of 20 km, a denser lower 



^Di, has been set equal to 18 km and 4 km for the Moon 
and Earth, respectively l llvanov et al.ll2boih . 
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crust down to 40 km, a peridotitic lithospheric 
mantle which is in turn characterized by an upper 
layer with an average density of 3.2 g/cm"^ (up to 
150 km), and a lower layer with an higher density 
(3.3 g/cm^). The impactor density has been set to 
2.7 g/cm'^. Regarding the strength, the value re- 
ported bv lAsphaug et al.l (jl996h for bulk siHcates, 
namely Yq = 2 x 10* dyne/cm^ was used both for 
the Moon and Earth crusts. However, we assume 
a linear increase of the strength, from zero at the 
surface up to Yq at the bottom of the heavily frac- 
tured layer for the Moon in order to take into ac- 
count the cohesionless regolith layer on the surface 
and the underlying megaregolith and anorthosites 
which are likely characterized by a progressive de- 
crease in fracture density with depth. A similar 
assumption was adopted also for the layer of sed- 
imentary rocks on the Earth. The strength at the 
Earth's surface was set to 0.3 xlO* dyne/cm^. 
Following our description of the Moon crust, its 
upper layers are mainly made of highly fractured 
materials, which behave like cohesive soils. There- 
fore, we set a sharp transition in the scaling law, 
from cohesive soil in the case of small impactors 
to hard rock for larger ones, at a projectile size 
of 1/20*'' of the thickness of the heavily fractured 
silicate layer (i.e. 0.5 km). For the Earth only the 
hard rock scaling law was used since only rocky 
layers have been assumed at the surface. 

3. The Model Production Function (MPF) 

We have now all the necessary inputs for com- 
puting the cumulative crater size distribution 
(hereafter the Model Production Function, MPF) 
for the target bodies. Let ^{D) be the crater 
differential distribution, we may write: 



A detailed comparison between t hese two pro- 



MPF(i:>) = / ^{b)db 



(9) 



where $ can be expressed in terms of and S 
(see appendix for details). In figure |4] our lunar 
model production function is shown for different 
scaling laws along with the Ncukum production 
function for the Moon (NPF). The NPF has been 
the traditional reference for dating purposes since 
the late 70s, with a few revisions in more recent 



times (jNeukum 1983; Neukum et al.ll200lD . An al- 
ternative production func tion, the so-calle d HPF , 
has also been proposed bv lHartmann et ah (1981). 



duct io n functions can be found in iNeukum et al.l 
(j200l[ ). For our purposes, we recall that the NPF 
shows an overall agreement with the HPF, even 
if some dis crepancies are presen t. Nevertheless, 



as stated in lNeukum et al.l (|200ll ). the NPF-based 
chronology represents the current best interpreta- 
tion for the lunar cratering chronology. Moreover, 
since the NPF extends to a wider size range than 
the HPF, the former is more suitable for a thor- 
ough comparison with the MPF. For these reasons, 
here we provide detailed comparison between the 
MPF and the NPF. In the following paragraph, 
we briefly recall how the NPF has been derived. 
There are several versions of t he NPF. In this pa- 
per w e refer to the old NPF ( Neukum fc Ivanov 



19941 ) and to the new NPF (jNeukum et al.ll200lh . 



In both cases, the underlying assumption is that 
the shape of the production function has been 
constant, with only the absolute number of craters 
changing over time. Therefore, the NPF was built 
by vertically re-scaling crater counts from areas 
with different ages, until they overlapped. The 
final curve was then expressed as a polynomial 
fit of the re-scaled data points. The strength 
of such procedure is that it is not model depen- 
dent because based on measured data. On the 
other hand this procedure, relying on the accu- 
racy of the measurements in multiple overlapping 
regions, is particularly prone to severe error prop- 
agation, especially at the large diameter end. In 
addition, in some regions the cumulative distri- 
butions of craters can be altered by surface pro- 
cesses like sporadic magma effusions, ejecta from 
other craters etc. The subsequent effects of su- 
perposition of geological units of different ages, 
normally not recognizable on remote sensing im- 
ages, can lead to mixed crater size distributions 
( Neukum fc Horn 19761 ) . Figure [4] clearly shows 



how the change of crater counting on single regions 
may alter the final shape of the NPF. In partic- 
ular, the discrepancies between the old and the 
new NPF are due to an adju stment of the crate r 
counts for Orientale Basin ( Ivanov et al. I2OOII) . 



A few considerations can be drawn from fig. [H 
The first one is that different scaling laws produce 
differences on the MPF as large as a factor of ten 
for large craters. We find that the MPF produced 
with the S&S scaling law is not in good agreement 
with the NPF throughout the whole size range. 
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When applying the same scahng law to t he ob - 
served NEO population, IStuart fc Binzell (|2004l ) 
found a better match to the new NPF for craters 
> 10 km. This discrepancy is mainly due to the 
different value of the intrinsic probability that 
Stuart and Binzel used in their model (namely, 
1.50 • 10~^ yr~^)- However, we believe the Bottke 
et al. determination of the intrinsic collision prob- 
ability to be more accurate since it relies on a bet- 
ter treatment of the gravitational focusing effects 
(A. Morbidelli, pers. communication). The main 
difference between the H&H and the I&S scaling 
laws takes place at crater sizes smaller than 10 km. 
This can be traced back to our choice of using the 
cohesive soil scaling law for small crater sizes in 
the H&H scaling law which causes an inflection 
point in the MPF. This transition in the target's 
crustal properties produces an S-shaped feature 
in the MPF, which is qualitatively similar to the 
one found in the NPF. This finding would suggest 
that the observed S-shaped feature in the NPF is 
caused by physical properties of the target, rather 
than reflecting the shape of the siz e distribution 
of the imp actors, as proposed by IWerner et al 



(|2002l ) and llvanov et all (|2002l ). Actually, a weak 
wavy feature in the relevant size range is present 
in the hn{d) which, however, is smeared out when 
passing to the MPF. However, in order to finally 
assess this issue, a better estimate of the observed 
impactor size distribution for sizes < 0.1 km is 
needed. It should also be mentioned that the de- 
tailed shape of the NPF somewhat depends on the 
fitting procedure used to construct the NPF itself, 
as shown by the mismatch between the old and 
new NPF for large craters. 

In conclusion, despite of all the caveats mentioned, 
it is a remarkable result that the present MPF ob- 
tained with the H&H scaling law and the NPF 
(both new and old) are similar within a factor of 
two throughout the whole size range considered, 
that is from 0.01 km to a few hundred km. On 
the other hand, the S&S and I&S scaling laws fail 
to accurately reproduce the absolute number of 
craters per km^ per yr, with a maximum devia- 
tion of a factor of 10. The good agreement of the 
absolute density of craters per yr derived from the 
observations (either the NPF or the HPF) and 
the MPF with the H&H scaling law is therefore a 
clear indication that the latter scaling law is the 
more appropriate in describing the lunar crater- 



ing. This agreement is particularly interesting, 
considering that the NPF and MPF are derived 
from completely independent methods: one being 
based on crater counting, and the other being the 
result of theoretical modeling. For these reasons, 
we decided to restrict the following analysis to the 
H&H scaling law. We refer to this MPF as to the 
nominal model. 



3.1. The reliability of the nominal MPF 

In this section we explore the effects that the 
various parameters involved in our model have 
on the nominal MPF described in the last sec- 
tion. In particular, we analyze the scaling law, 
the NEO size distribution and the assumed den- 
sity and strength profiles. 

The H&H scaling law depends on three parame- 
ters (namely k, v and [i) which have been derived 
by best fit to laboratory data. We applied a varia- 
tion of 10% around their nominal values, and com- 
puted the resulting MPF. It results that fc and v 
have a negligible effect on the MPF. Changing /i of 
±10%, the MPF shifts vertically of about a factor 
of ±2, respectively. 

A detailed comparison of the modeled and ob- 
served NEO size di stributions can be found in 
Bottke et al.l (|2005al lbl) . For our purposes, we no- 
tice that they basically overlap for d > 1.5 km. 
For smaller impactor sizes, the situation is less 
clear since different surveys produce slightly differ- 
ent observed distributions. Anyway, we can safely 
state that we have a maximum deviation between 
the observations and the model distributions of 
about a factor of 1.8 (being the observations lower 
than the model). 

The scaling law is weakly dependent on the den- 
sity, so a change in the numerical values affects 
the MPF in a negligible way. Notice also that 
craters larger than about 0.1-0.2 km form in grav- 
ity regime, hence the rock strength is not impor- 
tant at those sizes. For smaller sizes, the shape of 
the MPF may depend on the local strength. The 
density and strength profiles adopted in this work 
are meant to be average profiles for the Moon. 
Local variations to the mean values may slightly 
change the MPF at small crater sizes. 
A more important parameter is the assumed depth 
for the transition from cohesive to hard rock scal- 
ing law. A variation in this value causes a shift of 
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the inflection point in the MPF within the range 
1 < Z? < 10 km. 

4. A new chronology for the Moon 

The MPF can be used as a reference to develop 
the lunar chronology. For this purpose, we used 
literature crater counts from regions with known 
radiometric ages as calibration regions. A list of 
the regions used is reported in table [T] 
In order to obtain the lunar chronology, we first 
assumed that the MPF shape remained constant 
over time. For each calibration region we deter- 
mined the proportionality factor which gave the 
best fit to the data by minimizing the x^. Fi- 
nally, we determined the crater cumulative num- 
ber at 1 km (iVi). The fitting procedure used is 
particularly important, as it directly affects the 
chronology. For this reason, we performed several 
tests in order to assess its stability. A major as- 
pect consisted in choosing the weights given to the 
measured data points. It is common practice in 
crater counting to assign each data point an error 
bar corresponding to the square root of the num- 
ber of counts, under the assumption that the mea- 
surements follow a Poisson distribution. However, 
in the case in which significant systematic effects 
are present, this procedure leads to underestimate 
the total error. This problem is particularly evi- 
dent when dealing with small craters that are at 
the limit of the image resolution. Because of their 
large number, the statistical error for these craters 
is comparatively small. On the other hand, due to 
the difficulty of positively identifying craters at the 
limit of the resolution, these measurements could 
be affected by a significant bias. This effect may 
affect the A^i values up to a factor of 3-4 in some re- 
gions and consequently have a non negligible effect 
on the chronology. After several tests, we decided 
to use the following weights: 



af = N + kN^ 



(10) 



where k, was set to 0.5. In this way, the total er- 
ror {at) for points with very low statistical a were 
increased, while those with large statistical error 
(i.e. low number of craters) were basically left un- 
changed. 

In order to validate this procedure, we applied 
it to the old NPF, being able to reproduce the 



Neukurn fc Ivanovl (|l994r ) results within their er- 
ror bar^. 

4.1. Moon cratering 

In this section we present the results of the fits 
applied to the regions used to calibrate the lunar 
chronology (see tab. [T]). The fitting procedure 
was successful for most of the regions. In sev- 
eral cases a very good fit was achieved through 
the measured range of the cumulative distribu- 
tions, indicating an overall agreement between 
the shape of the MPF and the data. However, 
some discrepancies occurred for the highlands (for 
D > 100 km). Mare Crisum and Mare Tranquil- 
litatis. In the case of the highlands (see fig. [5]), 
the discrepancy seems to suggest a different shape 
of the impactor size distribution in the past, with 
a larger fraction of medium-to-large impactors 
(see section [5] for further details). In the case 
of Mare Crisum and Mare Tranquillitatis regional 
phenomena can be responsible for the observed de- 
viations. In particular, in both cases several con- 
secutive magma flows have partially covered the 
population of small crater s, affecting the observed 
cumulative distri bution i Neukum fc Horn Il976t 



Bovce et al.lll977[ ). Indeed the absolute radiomet- 



ric ages derived from All samples (Mare Tran- 
quillitatis), span from 3.5 to 3.85 Gyr, whereas 
the ages of the Luna 24 basalts (Mare Crisium) 
are clustered into at least three different peaks 
at 2.5, 3.3, 3.6 Gyr, respectivelv feirck fc Allegre 



1978[ IStoffler fc Rvdeil200ltlFernandes fc Burgess 



20051 ). In presence of such surface phenomena, the 
fit must be constrained on the basis of geological 
considerations, as the automatic fit may be mis- 
leading. In part i cular, as formerly suggested by 
Neukum fc HornI (|1976^ . the crater size- frequency 



distribution derived from these regions leads to a 
composite curve in which smaller craters reflect 
the younger radiometric ages, while the larger 
ones correspond to older ages. Thus, for Mare 
Tranquillitatis we forced the MPF to overlap the 
cumulative curve at 0.5 km for the age of 3.55 Gyr 
(referred to as young) ^ and at 0.9 km for the age 
of 3.72 Gyr (referred to as old). Similar consider- 
ations hold also for Mare Crisum where the MPF 



^Here we used the old NPF inst ead of the new NPF be cause 
even the latest chronology by iNeukum et al.l ll200lh uses 
the old NPF. 
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was fitted to craters with diameters in the range 
1 to 2 km (details on the fit of the caHbration 
regions are reported in the online material). 
Another effect influencing the accuracy of the 
crater statistics is the possible presence of sec- 
ondary craters. Even though authors publish- 
ing crater counts pay particular attention to the 
identification and rejection of secondary craters, 
in some cases confusion with pri mary craters i s 
difficult to avoid. According to llvano 
secondaries are likely negligible on young regions 
because they did not have the time to accumu- 
late in large numbers. On the other hand, older 
regions, namely the highlands and maria, may 
show secondar y craters caused by subsequent 
large impacts (IWilhelingI 19761 ). Nevertheless, 
Neukum fc IvanoX i 19941 ) concluded that the sec- 
ondaries are likely not relevant. For the present 
analysis, we assumed that the observed crater 
cumulative distributions are resulting from the 
primary flux of impactors. 

4.2. Earth cratering 

Cratering studies are much more difflcult on the 
Earth than on the Moon, due to strong resurfacing 
processes and consequent alteration of the crater 
morphologies. Nevertheless, Earth craters are im- 
portant because they can be studied in detail and 
precisely dated. Therefore they can be used to 
further constrain the flux of impactors, in partic- 
ular for young ages, that are scarcely represented 
in lunar data. In this work two data sets of terres- 
trial craters have been considered. The flrst one 
consists of all large craters {D > 20 km) found 
on the North A merican and Euroasiat ic cratons, 
as presented bv iGrieve fc Denc3 (|l979t ) (updated 
with few new findin gs from the Earth Im pact 
databas^. Following Grieve fc Pence! ( 19791 ). we 
have adopted for these two cratons an age of 
0.375 Gyr. Some sectors of both cratons were 
probably exposed to the meteoroid fiux already 
from the early Ordovician (0.45 Gyr), whereas 
other s only from 0.300 Gyr ago (jGrieve fc Pence 
19791 ). The choice of 0.375 Gyr is further sup- 
ported by the age of the older craters that, with 
only few exceptions, are around this value. 
The second data set represents a subset of the 
first one, including only craters with radiometric 



ages less then 0.120 G y r. Th is choice has been 
made following Grieve! ( 1993f l. in order to over- 
come biases due to erosion and burial of craters 
on the Earth's surface. Both cumulative distri- 
butions were fitted using our Earth MPF, and the 
obtained Ni was rescaled to the lunar case by con- 
sidering the gravitational focusing and the differ- 
ent relationships between impactor and crater size 
on the two bodies (see fig. [6]). The resulting val- 
ues for A^i are shown in tab. [1] with more details 
supplied in the online material. 

4.3. MPF chronology 



The lunar chronology has been developed on 
the basis of the relationship between the derived 
Ni values and the absolute ages. However, this 
task is not trivial since only 18 measurements are 
available, 12 of which are older than 3 Gyr. More- 
over, the region from 1 Gyr to 3 Gyr lacks com- 
pletely in data points. Therefore a very accurate 
chronology cannot be expected, in particular for 
ages younger than 3 Gyr. In this work the radio- 
metric ages proposed by IStoffler fc R,-^ (|200lh 
were used. 



In ord er to describe the lunar chronology, iNeukum 
(|l983l ) suggested the following approximation: 



iVi = a{e''* - 1) + ct. 



(11) 



This function assumes a linear relationship be- 
tween A^i and the age (t) for fei <C 1, while for 
bt ^ 1 Ni increases exponentially. Although this 
is a rather simple description of the data, it accu- 
rately fits the data points in the Neukum chronol- 
ogy. We fitted the same function to our MPF- 
based data points, with the results shown in fig. 
HI The derived coefficient are: a = 1.23 x 10"^^, 
b = 7.85, c = 1 . 30 X 10~^. For comparison, 
Neukum fc Ivanovl (|l994[ ) obtained: a = 5.44 x 
10-14^ 6 ^ 6.93, c = 8.38 x lO "'^, which are quite 
similar to our rcsull[£|. The Neukum fc Ivanov 



^As previously discussed, old regions with cumulative count- 
ing extending to sub-km sizes may be affected by sec- 
ondaries. These regions are: Descartes Formation, Mare 
Imbrium, Mare Fecunditatis, Taurus Littrow Mare, Mare 
Tranquillitatis. Therefore we also performed a best fit 
excluding these regions. The resulting coefficients are 



3.07 X 10" 



7.63, c 



1.33 X 10" 



and the 



*http://www.unb.ca/passc/ImpactD at abase/. 



resulting chronology curve is basically overimposed to the 
one obtained using all the regions. Therefore the secondary 
craters -if presents in the cumulative distributions- do not 
affect the chronology. 
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(|l994l ) chronology curve is also plotted in fig. [71 
While the derived best fit is in overall agreement 
with the nominal MPF-based data points, there 
are some deviations that mainly concern young 
data points, which lie systematically above the 
best fit curve. It is remarkable that in our chronol- 
ogy, the proposed equation fits well all the points 
older than Coper nicus (which was an ou tlier in 
the chronology of iNeukum fc Ivanov 1994I ). while 
it underestimates -by a factor of 2- all the points 
younger than Copernicus (see fig. [T]). 

4.4. Non-constant flux in recent times? 

As mentioned, the MPF-based chronology is 
reasonably in agreement with that derived from 
the NPF. The discrepancies observed for young 
regions could be due to uncertainties involved in 
our model. In particular, a concern regards the 
slope of the impactor cumulative size distribution 
for small dimensions which is poorly constrained 
by observations. Errors in the slope in this size 
range may result in inaccurate estimates of A'^i for 
the young lunar regions. The slope of the NEO 
model for impactor diameters smaller than 0.1 km 
(which roughly corresponds to 1 km crater size) is 
-2.6, which becomes a slope of -3.1 in the MPF. 
This slope is very similar to that of the NPF for 
the same size range, which is equal to -3.2. This 
slight difference in slope can account for half of 
the observed difference between our and Neukum's 
chronology. We underline that such slope varia- 
tions may affect the Cone, Tycho and North ray 
craters, but not the terrestrial crater counts since 
they have I? > 20 km. Therefore it is remark- 
able that terrestrial craters lie on the same linear 
trend of the young lunar data points (fig. [7l left 
panel). The NEO population model adopted here 
is in overall agree ment with the availab le observed 
data (see fig. 16 of lBottke et al.l[2005b( ). For small 
dimensions, however, the real NEOs size distri- 
bution is poorly known. The only observational 



const raints are the small fireballs (IHallidav et al 



19961 ) and bolides ([Brown et al.ll2002f ) detected in 
the E arth's atmosphere. A careful check at fig. 
16 of [Bottke et al.l ([2005bl ) shows that the NEO 
population model slightly overestimates -by a fac- 
tor of about 1.5- the predicted number of NEOs 
from bolide data. This occurs for the impactor size 
range responsible for the cratering observed for the 
Cone, North Ray and Tycho craters. Based on the 



data presently available, it is very difficult to un- 
derstand whether this difference is real or rather 
due to the involved uncertainties. Nevertheless, 
if we maintain the present number of NEOs at 
d = 1 km, we obtain that the Ni values corre- 
sponding to the younger lunar regions would be 
further increased. 

Other possible sources of uncertainties may arise 
from the scaling law. Nevertheless, variations in 
the parameters produce a vertical shift of the MPF 
and not a change of the shape, therefore this shift 
alone cannot affect the determination of iVi values. 
The only way to obtain a change in the Ni values 
with respect to the nominal model, is to introduce 
a slope variation in the MPF for D < 1 km. This 
would be the case, if the inflection point would 
slide below D ^ 1 km. This would imply a very 
low thickness for the fractured layers which seems 
improbable, according to our understanding of the 
lunar crust. Nevertheless, if this would be the 
case, Ni for the young lunar regions would be fur- 
ther increased. 

It is also possible that the MPF-based A^i values 
are accurate , and th e simple relationship propo sed 
by [NeukunJ ([l983[ ): [Neukum fc Ivanovl ([l994h is 
not adequate to describe the data. This would be 
the case, for example, if the impactor flux had not 
been constant during the past ~3 Gyr. Due to the 
lack of data points in the range 1-3 Gyr, we are not 
able to study in detail possible variations in the 
flux. However, fig. [7|seems to suggest that a recent 
phase (age < 0.4 Gyr) characterized by a constant 
impactor flux was preceded by a phase (between 
0.8 and 3 Gyr) with a lower and nearly constant 
impact rate. A single event placed around 0.4 Gyr 
ago and lasting for 0.2-0.3 Gyr that would increase 
the impactor flux by a factor 2-3 would explain 
the observed data. Such a recent change in im- 
pactor flux would also have affected the older re- 
gions, but here the relative increase in the number 
of craters would be very low compared to the accu- 
mulation of craters over 3 Gyr or more, hence this 
effect would be negligible. This scenario is quali- 
tatively in agreement with the recent suggestion of 
an increase in the number of impactors in the in- 
ner Solar System due to the formation of dynam- 
ical families -as a results of cata strophic disrup- 



tion of asteroids- in the main belt ([Nesvornv et al 



2007[ : [Bottke et al.[[2007l) . In particular, the Bap- 



tistina and Flora families are estimated to have 
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formed about 140 Myr and 500 Myr ago, respec- 
tively. Their proximity to the resonances led to 
an increase in the member of the NEO population, 
starting few tens of Myr after their formation, as a 
consequence of the Yarkowsky-effect-driven decay 
of the semimajor axis. In the case of Baptistina 
the maximum flux was reached about 60 Myr af- 
ter its formation (i.e. 80 Myr ago). Therefore this 
flux may have affected the Cone, North Ray and 
Tycho craters counting, and also the young ter- 
restrial craters. 

Concerning the formation of Flora, similar consid- 
erations may hold. In addition, the Flora family 
has been connected to the spike in the infall of L- 
chon drite meteorit es that occurred about 470 Myr 
ago ( Bogard 19951 ). A similar spike has also been 
detected in the lunar glass spherul es and -to a 



lesser extent- in lunar impact clasts ([Culler et al 



20001: ICohen et al.ll2005^ . These spikes may have 



affected mostly the terrestrial cratering and the 
Copernicus counting. Although such enhance- 
ment of the fl ux was probably form ed by sub- 
km impactors (jHartmann et al.l 120071 ). the spike 
recorded in the lunar melt clasts implies also the 
existence of a number of km-sized impactors. This 
enhancement in the km-sized impactors may ex- 
plain why the Earth cratering for _D > 20 km 
seems not to be depleted b y erosion proc esses (fig. 
[71 left panel), as expected ( Grieve! 1993 ). 



4.5. On the early flux 

One of the major -and mostly debated- open 
issues is related to the cratering rate in the early 
times after the formation of the Moon. Argu- 
ments have been proposed in favor of a rapid 
and smooth decrease in the impact rate from 
the formation of the Moo n till about 3.5 Gyr 



ago (jHartmann et al.l 11981 



, 1983), fol- 

lowed by a constant imp act rate (fNeukum 1983t 
Neukum fc Ivanovl Il994[) . On the other hand, 
studies of impact melts led some researchers (e.g . 



Ryder 



01 impad 
19901: ISt5 



Stoffler fc Rvded 120011 : ICohen et al 



20001) to support the idea of an intense lunar bom 



bardment about 3.9 Gyr ago which lasted some 
100-200 Myr. From a dynamical point of view, 
such short phase of intense bombardment has been 



^See the discussion in lBottke et al.l | |2007|) about the prob- 
abihty of Tycho itself being formed by an impact with a 
Baptistina fragment. 



recently explained in the more general context of 
the e arly stages of the e volution of the Solar Sys- 



tem ( Gomes et al. 20051 ). This scenario has also 



found some observational evidence in the work of 
( 20051) . 



Strom et al 



Despite these works, the LHB hypothesis has not 
found an unanimous co nsensus yet. For instance. 



Hartmann et al.l (|2007l ) raised doubts about the 



inter pretation of lunar melts and glasses data; 
while lNeukum fc Ivanovl (|l994l ) argued against the 
LHB on the basis of the smooth behavior of their 
lunar chronology curve. 

If the intense LHB did take place, it should have 
left some traces in the chronology curve, which 
therefore could be used to obtain constraints on 
the early stages of the cratering. Unfortunately, 
there are only 5 measured regions having ages 
older than 3.85 Gyr, and therefore it is not easy to 
draw firm conclusions. One of the most important 
region s is the Nectaris Basin. iNeukum fc Ivanov 
(|l994l ) used a radiometric age of 4.1 Gry. The re- 
sulting best fit for the lunar chronology curve is 
therefore very close to all the old regions (within 
errors) and has a smooth behavior in early times, 
reflecting a possible smooth decay in the impact 
rate (see fig. [7]). However, using the new esti- 
mate of 3.92 Gyr for N ectaris Basin proposed by 
IStofHer fc B^Me^ (|200lh this point moves consider- 
ably far from the best fit (this is the age we used 
in our chronology) . This results in a change of the 
slope of the chronology curve, which in turn is in 
favor of the LHB. We underline that in the present 
analysis we used the present NEO size distribu- 
tion. However, it is very likely that at the time 
of the LHB the impactor size distribution had a 
different shape. At the present stage of knowledge 
a firm conclusion from the chronology curve can- 
not be obtained. Nevertheless, we will add few 
considerations on this important point in the next 
section. 

5. A time-dependent MPF 

In this section we develop a time-dependent 
MPF. In doing so, we clearly show the versatility 
of our approach, as opposed to the conventional 
methods (e.g. the NPF), where the production 
function is generally assumed to be constant over 
time. Modeling a time-dependent MPF is quite 
complex, and in this section we will only pose the 
basis of the problem, deferring thorough investiga- 
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tions to a further analysis. In general terms, there 
are two possible reasons for such time-dependent 
behavior: The first one is related to a non sta- 
tionary flu50j (t'id-TVjt); the second one to crater 
erasing processes. 

A motivation for studying the non stationary flux 
is provided by the dynamical processes responsi- 
ble for the flux. In early times (just after the for- 
mation of the terrestrial planets and the Moon, 
some 4.5 Gyr ago) the inner Solar System was 
still populated by the leftover planetesimals, which 
were rapidly (10-100 Myr) cleared-up by the latest 
stages of the accretion process. A second source 
of impactors was represented by the highly per- 
turbed bodies (e.g. through the sweeping reso- 
nance mechanism) due to the formation and mi- 
gration of Jupiter. At the present time, the flux 
is mainly sustained by the slow decay of MBAs 
into the resonances. All these three stages have 
their own characteristic 4){d,v). For sake of sim- 
plicity, here we focus on the last two mechanisms. 
The second process is size dependent, while the 
first one is not, therefore the shape of the im- 
pactor size distribution changes over time. On the 
other hand, the impactor velocity distribution is 
expected to remain approximately constant over 
time, as in both cases the transport mechanism 
to the NEO space is the same. Therefore we may 
write 4>{d,v,t) = h{d,t)f{d,v). 
The crater erasing on large bodies depends mainly 
on two major effects: the superposition of craters 



and the local jolting (jO'Brien et al.ll2006f) . The 



first one basically accounts for the overlapping of 
craters as the surface crater density increases. The 
second one considers the erasing of craters due to 
local seismic jolting triggered by the formation of 
large craters. Both effects can be modeled by con- 
sidering a vari ation in the number of craters in a 
given size bin ( O'Brien et al. 2006h . Let £ be the 



ratio between the actual number of craters consid- 
ering the erasing and the total number neglecting 
the erasing, we can write: 



MPF(D,i) 



^{D,t)£{D,t)dD (12) 



D 



^We stress once again that in our approach what is relevant 
is the shape of the impactor size and velocity distributions, 
since the absolute calibration of the model is done by the 
reference lunar regions. 



where setting t — and £ — 1 reduces to 
the MPF used in previous sections (see appendix 
for further details). Here we limit ourselves to 
considering these effects separately. Let us con- 
sider first the crater erasing mechanism, assuming 
(j) = (j)(d,v,0) as in previous sections. The re- 
sults are shown in fig. [51 It is clear that as the 
age increases the MPF becomes flatter for small 
dimensions. This is the well-known saturation 
process of heavily cratered surfaces. This affects 
the derived Ni values. To study possible varia- 
tions in the rate of impacts over time, we adopted 
the following procedure: first the cumulative dis- 
tributions for the calibration regions were fitted 
with the MPF corrected for the erasing, then the 
obtained Ni was rescaled to the total number of 
craters that occurred in the region (i.e. the actual 
number of impactors). In fig. [9] the effects of the 
erasing on the chronology are shown. Although 
the erasing considerably changes the MPF for the 
old ages (see fig. |5]), its effects on the chronology 
are not important. 

Let us now deal with changes in the impactor size 
distribution and no erasing {£ = 1). Following 
the previous discussion, we assumed that before 
the end of the LHB (~ 3.8 Gyr ago) the size dis- 
tribution h{d) resembled closely that of the main 
beliEl (hm), while afterwards it becomes similar 
to the present one, that is This assumption 
is also encouraged by the observation that hm 
provides a much better fit to the crater distribu- 
tion of the highlands and Nectaris Basin than /i„ 
does (see online material). With respect to the 
previous analysis (see fig. [7]), the only difference 
regards two data points, namely the highlands 
and Nectaris Basin. For both regions we obtain 
a lower Ni value, with the reduction being more 
pronounced for the highlands. This results in a 
net change in the slope of the chronology in the 
two oldest points, suggesting a change in the rate 
of impactors. Therefore, if the MBA size distri- 
bution is used for earlier times on the basis of the 
improved fit with both the highlands and Nectaris 
Basin, the chronology plot would suggest the ex- 
istence of the LHB. These conclusions are weakly 
affected by the uncertainties in the model (size 
distribution, scaling law etc.) because the MPF is 



'For sake of simplicity, we used the present MBA size dis- 
tribu tion although it may have had a different shape in the 
past llBottke et al.ll2005al lbh. 
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well defined for the large crater sizes relevant in 
this context. 



6. Mercury's chronology 

The impactor size and velocity distributions for 
Mercury have been computed in a similar way as 
done for the Moon (figs [TOl fTT|) . We assumed den- 
sity and strength profiles similar to those adopted 
for the Moon. The resulting MPF is shown in 
fig- [m (right panel), along with the overp lot of 
the NPF for Mercury ijNeukum et al.ll2001bh . The 



two production functions nicely overlap at the ex- 
tremes, that is for D < 0.3 km and D > 40 km. In- 
side this range, however, the two curves are quite 
different, with a maximum deviation of nearly a 
factor of 5 at - 3 km (fig. [H left panel). The 
exact behavior of the MPF in this range depends 
on the assumptions on the transition in the scaling 
law between cohesive soil and rock (see discussion 
in sect. 12. 2p . Nevertheless, whatever reasonable 
choice is made for the scaling law parameters, our 
MPF is remarkably different from the NPF in this 
range. 

As an example, we report here our age determi- 
nation of the Chekhov basin. Crater counting 
was performed on the total ejecta blanket, over 
an area ex ceeding 4 ■ 10^ km^ . A ccording to the 
analysis of Neukum et al. ( 2001bl ). Chekhov basin 
has a derived age of 4.05 Gyr, which makes it 
one of the oldest regions on Mercury. In fig. [T3] 
we report our MPF-based age estimate, obtained 
considering different scenarios, namely using the 
NEO and MBA size distributions. In the nominal 
case (i.e. MPF=MPF(i:), 0); £ = 1), we obtain 
an age of 3.94 Gyr, which is s l ightly y ounger than 
that found bv lNeukum et al.l (|2001bl ). The corre- 
sponding of the best fit is 0.8. Note that, in 
comparison to the NPF (x^=1.5), the MPF more 
closely fits the observed crater cumulative distri- 
bution. When considering also the crater erasing, 
the age becomes 4.05 Gyr (x^=0.6), and the shape 
of the MPF is basically unchanged. In the case 
of the MBA size distribution and no erasing, we 
obtain 3.97 Gyr (x^=0.6). Notice that in the case 
of MBA, if we include crater erasing the MPF 
reaches the saturation prior to achieving the best 
fit, therefore it is not possible to derive a precise 
age estimate. In conclusion, the MPF shows a 
better fit to the data points than the NPF. How- 



ever, due to the measurement errors and to the 
limited crater size range observed, it is not pos- 
sible to conclude whether the NEO or the MBA 
size distribution is more suitable for fitting the 
Chekhov basin. 



7. Conclusions 

In this paper a new model for the chronology 
of planetary surfaces in the inner Solar System 
has been introduced and applied to the Moon and 
Mercury. Concerning the Moon, the main findings 
are the following: 

• The nominal MPF differs from the NPF for 
less than a factor of 2 throughout the ob- 
served range of dimensions. A similar degree 
of agreement is obtained for the chronology; 

• Despite the good agreement between MPF 
and NPF, there is a systematic misfit of 
the assumed linear branch of the chronology 
function with the data points for ages more 
recent than about 0.4 Gyr. It is remark- 
able that the Earth craters (which have D > 
20 km) are aligned with the Tycho, Cone 
and North Ray craters (having D < 0.4 km). 
Although the misfit may be partially due to 
the uncertainties present in the model, other 
causes, such as for instance a non-constant 
flux in recent times due to the formation 
of MB families (like Baptistina and Flora) 
might be possible; 

• The MPF was computed for two differ- 
ent impactor size distributions, namely the 
present NEO and MBA populations. In 
this way we were able to study whether the 
cratering production function changed over 
time. Both the NEO and MBA-based MPFs 
are able to fit the cumulative distributions 
having ages < 3.9 Gyr. This is mainly due 
to the fact that such regions have craters 
smaller than ~ 10 km, where the MPFs are 
very similar. On the other hand, the old- 
est regions (highlands 4.35 Gyr; Nectaris 
Basin 3.92 Gyr) are best fitted using the 
MBA size distribution. This suggests that 
the impactor size distribution around 4 Gyr 
ago resembled the present MBA distribu- 
tion, while today it is closer to the present 
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NEO size distribution. Alternatively, it is 
also possible that comets had in the past a 
non negligible contribution. This point de- 
serves a specific analysis and will be subject 
of a future work; 



Concerning Mercury, the MPF differs from the 
NPF in the size range from 0.3 km to 40 km. The 
maximum difference -of about a factor of 5- oc- 
curs at ^ 3 km. This discrepancy is therefore 
potentially important for dating purposes, in par- 
ticular when measurements for small craters will 
become available (e.g. from MESSENGER and 
BepiColombo). Outside this range, the NPF and 
MPF are very similar. We show the results of our 
model for the case of the Chekhov basin. In the 
nominal case (MPF(£', 0), 5 = 1) we obtained the 
ag e of 3.95 Gyr, which c ompares to the 4.05 Gyr 
of lNeukum et al.l (|2001bl) . 
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A. The Model Production Function 



The cumulative distribution of craters (the so-called production function) can be evaluated starting from 
the differential distribution of the incoming flux, namely (f){d, v, t). Here we consider the general problem of a 
non stationary flux, therefore the time dependence of the function 0. Each impact produces a crater, with a 
size specified by the scaling law, D = S{d, v,p). Let be the differential distribution of craters, namely 

the number of craters per unit of crater size per unit surface. Such distribution is obtained by integrating the 
differential distribution of impactors over the one-dimensional domain specified by the scaling law, namely 
7D : S{d,v,p) = D. Let ct be a parameter along the curve jjy, we finally have: 



da 



da 



(Al) 



In the presence of erasing processes, the actual number of measured craters in any given size bin dD is 
not equal to the number of impactors. Let £{D,t) be the ratio of the final number of crater erasing included, 
with the total number (i.e. erasing excluded), we finally obtain: 



MVY{D,t) 



^{D,t)£{D,t) dD 



(A2) 



Following the work of lO'Brien et al. I (l2006l) . the function £ can be expressed as 



£{D,t) = 



A 



(A3) 



where A is the total effective area of craters accumulated per unit area. Therefore A depends on the 
age since older regions have a large number of craters accumulated and hence a larger A with respect to 
younger regions. A also depends on crater diameter since small craters are more frequent than larger ones. 
As a results, the counting of small craters on old regions are strongly affected by the erasing (see fig. [8]). In 
calculating A, we considered the superposition of craters and the local jolting, and neglected global jolting 
and cumulative effects of seismic shaking, since they are not relevant for large bodie s like the terr estrial 
planets. The parameters involved in the estimate of A have been derived from O'Brien et al. (|2006l ). The 
implementation of the erasing processes in our code has been tested accurately using asteroid crater counting 
(Gaspra, Ida, etc). Our results agree with those published. Nevertheless, the purpose of the present analysis 
is simply to show the first-order effects of the erasing, and a detailed evaluation of the erasing on planetary- 
sized bodies is left for further analysis. 
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Table 1 

Description of the lunar and terrestrial calibration regions used in this work. 



Region I Ni (NEO)t Ni (MBA)+ Age* 







(km 








(km 






(Gyr) 


Highlands'^ 


7 


851 


10" 




2 


018 


W 




4.35 


Nectaris Basin 


1 


327 


10" 




6 


648 


10^ 




3.92 


lZ)cs cartes Format ion ^ ' ^ 


2 


490 


10" 


2 


2 


509 


10' 


-2 


3.92 


Imbrium Apennines^ 


1 


968 


10" 


2 


1 


931 


10" 


-2 


3.85 


Fra Mauro Formation 


2 


595 


10" 


2 


2 


672 


10' 


-2 


3.85 


Mare Tranquillitatis (old)'^'^ 


1 


836 


10" 


2 


1 


832 


10' 


-2 


3.80 


Taurus Littrow Marc^ 


1 


579 


10" 


2 


1 


585 


10' 


-2 


3.70 


Mare Tranquillitatis (young)^-^ 


9 


300 


10" 


3 


9 


357 


10' 


-3 


3.58 


Marc Fccunditatis^ 


3 


234 


10" 


3 


3 


257 


10' 


-3 


3.41 


Mare Imbrium^ 


5 


468 


10" 


3 


5 


526 


10' 


-3 


3.30 


Mare Crisium 


2 


335 


10" 


3 


2 


377 


10' 


-3 


3.22 


Oceanus Procellarum 


3 


683 


10" 


3 


3 


695 


10' 


-3 


3.15 


Copernicus Crater (cont. ejecta)^ 


1 


321 


10" 


3 


1 


337 


10' 


-3 


0.80 


Copernicus Crater (crater floor) *^ 


1 


348 


10" 


3 


1 


343 


10' 


-3 


0.80 


Terrestrial Phanerozoic craters^ 


1 


267 


10" 


3 


7 


655 


10' 


-4 


0.375 


Terrestrial Phanerozoic craters (young) ^ 


3 


835 


10" 


4 


2 


195 


10' 


-4 


0.120 


Tycho crater (cont. ejecta)'^ 


3 


391 


10" 


4 


3 


401 


10' 


-4 


0.109 


Tycho crater (cont. ejecta)* 


1 


644 


10" 


4 


1 


712 


10' 


-4 


0.109 


North Ray crater^ 


1 


389 


10" 


4 


1 


421 


10' 


-4 


0.053 


Cone Crater^ 


6 


970 


10" 


5 


7 


131 


10' 


-5 


0.025 



Note. — Where not explicitly quoted, the data are from INeukumI (I1983D and refer- 
ences therein. The plots showing the fits of the calibration regions are reported in the 
online material. 

^The derived Ni values both using the NEO and MBA size distributions. 

^Radiometric ages for the lunar regions are from table VI oflStoffler & Rvderl ( f200l|) . 
The ages of the terrestrial craters arc from [Grieve & Pence! 1 119791) ■ 

^All craters and basins. 

"^The counting f or Descartes For mation is limited to dimensions smaller than 1.2 km, 
while that used in INeukumI < |1983t) extends to larger sizes. 

'^The crater counting for Imbrium Apennines ( used in th i s wor k) is a sub sample of 
the crater counting from Imbrium Basin, used by INeukumI U983D . 

*^Original data from lGreeley & GauT^ ( 119701) . 

°For the Copernicus crater we have two distinct measure ments correspon ding to the 
crater floor and the continuous ejecta blanket. In fig. 31 of |Neukum| jl983l) . these sets 
were probably merged into a single distribution. We decided to keep these two sets 
separate, and to provide individual fits (see fig. [7]and the online material). 

^Te rrestrial Phanerozoic craters are derived from IGrieve &z Denc^ 

I 19791), and update d according to the Earth Impact Database 

i http://www.unb.ca/passc/ImpactDatabase/ ) . Only large craters {D > 20 km) 
on the North American and Euroasian cratons have been considered. For sake of 
completeness we report the list of craters used. In the North American craton: 
Beaverhead, Carswell, Charlevoix, Clearwater East, Clearwater West, Haughton, 
Manicouagan, Manson, Mistastin, Presquile, Saint Martin, Slate Islands, Steen River, 
Sudbury. In the Eurasian craton: Boltysh, Kamcnsk, Kcurussclka, Lappajarvi, 
Puchezh-Katunki, Siljan. 

^Only large {D > 20 km) and young (< 120 Myr) terrestrial craters. They are: 
Carswell, Haughton, Manson, Mistastin, Steen River, Boltysh, Kamensk, Lappajarvi. 

^'^Two distinct measurements were available for the Tycho crater. Both correspond 
to the continuous ejecta: in (h ) sm al l crate rs were counted, while in (i) large craters 
were counted. In fig. 31 of Ncukuinj ]1983l) these files were merged into a single dis- 
tribution. We decided to keep these two sets separate, and to provide individual fits 
(see fig. [TJ. Notice that [Neukuml | |1983l) , also reports 4 measured points in the range 
0.5-1.1 km that were not available to us. 

^Original data from lMoore et al.| ( 119801) . 
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^ These arc relatively old regions having eratcr counting extending to small sizes 
(< 1 km) and thus potentially affected by secondary craters. In order to check whether 
these regions affect the final MPF-based chronology, we also performed a best fit ex- 
cluding these points. The results of the fit is basically unchanged (see text for further 
details.) 
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Fig. 1. — Distributions of impact velocities on the Moon and on the Earth, for two different impactor sizes. 
The smallest (S) and largest (L) impactor sizes used in these simulations are 0.1 m and 72 km, respectively. 
The average impact velocities are 18. 6 and 20.0 km/s, respectiv ely. Notice that these average impact velocities 
are slightly lower then those used bv lNeukum fc Ivanov (1994), after applying the correction for the average 
impacting angle of 7r/4. The distributions are quite spread and therefore a large number of impacts occur 
at velocities considerably different from the average, thus affecting the final crater distribution (see text). 
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Fig. 2. — Model cumulative size distributions of impactors on the Mo on and on the Earth (after iBottke et al 



2002 ). An arbitrarily rescaled MBA cumulative distribution (after [Bottke et aL 2005a[ ) is also shown for 
comparison. The main belt distribution shows an S-shaped feature in the 0.2-2 km size range. Also the 
NEO distribution shows a qualitatively similar feature. The two distributions, however, have quantitatively 
different shapes, due to selection mechanisms taking place during the dynamical transport from the main 
belt to the near-Earth space. 
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Fig. 3. — Assumed density profiles for the Moon and the Earth. Left panel: lunar and terrestrial density 
profiles vs depth. Right panel: average density vs impactor size. The average density has been obtained 
by averaging the density up to a depth of 10 x the impactor radius. A similar consideration holds for the 
strength (see text for details). 




Fig. 4. — Comparison of the MPF obtained with different scaling laws, indicated by H&H, l&S, S&S (see 
section for details). The old and new NPFs are also shown. 
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Crater diameter (km) 

Fig. 5. — The plot shows the best fit of the crater-size cumulative distribution for the lunar highlands with 
both the MPF and the NPF. 




Crater diameter (km) Crater diameter (km) 



Fig. 6. — Detailed comparison of the lunar and terrestrial MPFs (using the H&H scaling law). 
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Fi g. 7. — The plot shows the Ni values for the calibration regions obtained using the MPF and those obtained 
by I Neukum fc Ivanovl ( 1994 ). The solid curves, namely the chronology curves, represent the best fit of the 
two sets of data points. The shadowed regions encompass a factor of ±2 around the MPF chronology curve. 
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O.I 1 10 

Crater diameter (km) 



Fig. 8. — Moon isochrones obtained with MPF(D,0), and including the crater-erasing (CE) mechanism 
The effect of longitude variations is also shown by the shadowed area for the 4.5 Gyr curve (wit h no erasing) 
The e xpected longitudinal asymmetry in the impactor flux has been modeled using eq. 1 of iMorota et al 



(|2005l ). The impact of such analysis on the A^i values is negligible, therefore the chronology is not affected. 
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Fig. 9. — In this figure the effects of using the size distribution of MBAs and the crater-erasing (CE) 
mechanism are shown for the old regions alone, since the young ones are not affected. The vertical strip 
corresponds approximately the LHB event. Notice t hat for Nectaris basin th e age from Stoffier &: Rvdeil 
(j200l[ ) adopted in this work differs from the one in iNeukum &: Ivanovl (jl994l ) . Moreover, the age of the 
highlands is k nown with large uncertain ty and may vary in the range 4.2-4.4 Gyr. Here we adopted the 
values used bv lNeukum fc Ivanov ( 1994 ). 
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Fig. 10. — The plot shows the model cumulative distribution of impactors for Mercury, in comparison with 
that of the Earth-Moon system. The MBA size distribution is also shown. 




Impact velocity (km/s) 



Fig. 11. — Impactor velocity distribution for Mercury, in comparison with the Earth-Moon system (see also 
fig. ©. 
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Fig. 12. — MPF for Mercury, in comparison with the NPF. The two production functions are in disagreement 
in the range from about 0.1 km to 40 km. 




Fig. 13. — MPF best fit of Chekhov basin and derived age. We also present the best fit using the MBA size 
distribution and crater erasing (see section [5|). 
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